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Abstract 

Background: A response surface analysis was performed to study the effect of the composition and feeding 
thermal conditions of ternary mixtures on the number of theoretical stages and the energy consumption of Petlyuk 
columns. A modification of the pre-design algorithm was necessary for this purpose. 

Results: The modified algorithm provided feasible results in 100% of the studied cases, compared with only 8.89% 
for the current algorithm. The proposed algorithm allowed us to attain the desired separations, despite the type 
of mixture and the operating conditions in the feed stream, something that was not possible with the traditional 
pre-design method. The results showed that the type of mixture had great influence on the number of stages 
and on energy consumption. A higher number of stages and a lower consumption of energy were attained with 
mixtures rich in the light component, while higher energy consumption occurred when the mixture was rich in 
the heavy component. 

Conclusions: The proposed strategy expands the search of an optimal design of Petlyuk columns within a feasible 
region, which allow us to find a feasible design that meets output specifications and low thermal loads. 
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Background 

Distillation, one of the unit operations employed for 
separation of vapour-liquid mixtures, is a very energy 
demanding process. The use of non-conventional distil- 
lation columns, such as those suggested by Petlyuk et al. 
[1], can generate significant savings in energy consump- 
tion and thereby significantly reduce the capital costs 
compared to conventional sequences. The use of ther- 
mally coupled columns has decreased the costs of equip- 
ment and energy consumption by up to 30% [2-8]. 

Figure 1 shows examples of separation sequences by 
distillation, both conventional (a and b) and non- 
conventional (c, d, e and f), for the separation of ternary 
mixtures. Configurations c) and d) are employed fre- 
quently in the oil industry for the separation of air by 
cryogenic distillation (for example, a side rectifier is used 
for the recovery of argon). The e) configuration is a fully 
thermally coupled column, known as a Petlyuk column 
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[4,5,7-11]. The f) configuration corresponds to the di- 
vided wall column [10], also known as Kaibel column; 
these columns are currently used in the industry [11]. 

The proposed methods for the design of thermally 
coupled columns involve a pre-design stage that consists 
in breaking the system into a set of equivalent columns. 
The results from the pre-design stage are used to deter- 
mine the initial design configuration of the Petlyuk 
column, which must be optimized in order to attain 
minimum energy consumption. 

Other works about the design of Petlyuk-type columns 
have reported case studies of hydrocarbon mixtures with 
an equimolar composition [5-8] or with a high concen- 
tration of some of the components [12-15]. However, the 
pre-design method reported fails many times (no conver- 
gence to a solution) or yields infeasible solutions (poor 
quality output streams) when applied to other mixtures or 
other feed conditions such as non-equimolar mixtures 
and mixtures in which no single component is more 
highly represented than others. These facts are demon- 
strated in this work. 

This paper proposes modifying the algorithm used in 
the pre-design of Petlyuk-type columns for distillation of 
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Figure 1 Conventional (a, b) and non-conventional (c, d, e and f) sequences for separating ternary mixtures. 



ternary mixtures so that the pre-design stage is always 
successful and guarantees feasible results. Separation 
fractions are established in the equivalent configuration 
by minimizing the total sum of the theoretical stages 
and restricting the composition of the output streams to 
a desired value. 

Methods 

The method used for designing Petlyuk-type columns 
consists of two stages: a) pre-design and b) optimal design. 

Pre-design 

This stage consists in creating a distillation sequence 
equivalent to the Petlyuk column. The equivalent config- 
uration for separating three components is shown in 
Figure 2; it contains a sequence of three columns. Each 
column of this sequence is designed using the Fenske- 
Underwood-Gilliland (FUG) procedure for multicompo- 
nent distillation with the shortcut method [4,16]. 

The FUG design method requires knowing the following 
about each column; a) feeding conditions, (established by 
the problem to be solved and by the interconnection be- 
tween columns); b) the operation reflux ratio (fixed to 1.2 
times the minimum reflux); c) pressure and type of con- 
denser (fixed according to the FUG method); d) recovery 
fractions of key components. Conditions (a) to (c) are 
established in the traditional manner of the FUG method 
as described by Henley and Seader [16]. However, the 
manner in which condition (d) is established is where the 
difference between the traditional pre-design method and 
the method proposed in this article lies. 

In all cases reported in the literature [2,3,13,15,17], the 
recovery fractions for all columns are set directly, with 




Figure 2 Equivalent configuration of a Petlyuk-type column for 
separating a ternary mixture. 
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Figure 3 Algorithms of the solution for pre-designing Petlyuk-type columns: a) Traditional method with a fixed fraction; b) method 
proposed in this work. 



values ranging from 95% up to 98% in both key compo- 
nents. In the method proposed in this work, the recov- 
ery fractions are established by an optimization problem 
in which the total number of theoretical stages of the 
equivalent configuration is minimized and the compos- 
ition of the output streams are restricted in order to 
reach their desired value. The target function of the 
optimization problem, which must be minimized, is pre- 
sented in Equation 1: 



Si (fu<,l J HK.l JlK.2 JhK2 -JlK3 -JhK^ 
3 



(1) 



where / is the recovery fraction of key light components 
(sub-index LK), or key heavy ones (HK), in the respective 
column (sub-indexes 1, 2 or 3); A/j is the number of the- 
oretical stages of the column i (1, 2 or 3) and N is the 
total number of theoretical stages of the complete sys- 
tem. The problem, set forth with only the objective func- 
tion, contains 6 degrees of freedom; however, 4 of them 
are equality restrictions that restrict the composition of 
the desired product. In this sense, the composition of 
the light component (component A) in the top product 
stream of column 2 is established at a desired value of pur- 
ity; the composition of the intermediate component (com- 
ponent B) in the bottom product stream of column 2, as 



To use the pre-design of Petlyuk column. 



Optimization. 



Objetive: to obtain the reflux ratio so that SSE is minimized 
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Restrictions of equality. 

Compositions of the output flows are similar to those wanted. 



Restrictions of inequality. R op >R m 



Optimization. 



Objetive: to determine the interconection flows that minimize 
the thermal load in the main column. 
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V inf > 0.0 



Figure 4 Algorithm for optimal design of Petlyuk-type columns. 
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well as the top product stream of column 3, are estab- 
lished at the same desired value of purity, while the 
composition of the heavy component (component C) in 
the bottom product stream of column 3 is established at 
a desired value of purity. In this work, a purity of 98% 
was established for each of the products; however, any 
other value may be applied and the method will still 
work correctly. This composition allows the streams 
produced to present a high level of purity. With these 
equality restrictions, the problem of optimization is left 
with only 2 degrees of freedom and the feasibility of the 
solutions is ensured. Inequality restrictions are estab- 
lished in order to ensure that the problem can be 
solved. In any column, the minimum reflux ratio must 
be greater than zero, the top recovery fraction of the 
light key component must be greater than the top 

Table 1 Whole experimental plan with 90 study cases 



recovery fraction of the heavy key component, and must 
have at least 3 stages. Figure 3 shows the algorithms for 
the two cases discussed. 

It has been found that in the traditional method, in which 
the separation fractions are specified, the composition of 
the output streams differs from the desired specifications, 
mainly when the feed composition is not equimolar. In 
contrast, the method proposed here ensures that the com- 
position of the output streams is equal to the one specified 
by the user, even though the composition and thermal con- 
ditions of the feed streams vary. 

Optimal design 

The information obtained from the equivalent configur- 
ation (theoretical stages, molar reflux ratio, interconnec- 
tion stream flows, flow of each component in each output 
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stream and the position of the feed stage) is used for mak- 
ing the initial estimate in the second design stage of the 
Petlyuk-type column. At this point, the most valuable in- 
formation of the proposed pre-design method is the infor- 
mation of the flows of the interconnection streams. 

Using rigorous methods, it is possible to determine at 
this stage the operating conditions that are necessary for 
the Petlyuk-type column to consume a minimum amount 
of energy. In our work, this stage performed in the same 
way as reported in the literature [1,7]. The algorithm is 
represented in Figure 4. The optimal design is developed 
by performing two optimizations: the reflux ratio of the 
main column is determined by the first optimization, and 
the flows of the interconnection streams, which minimize 



the thermal load of the Petlyuk column are determined by 
the second one. 

The value of the reflux ratio is obtained by minimizing 
the objective function presented in Equation 2, which 
represents the sum of the squares of the errors between 
the calculated and the desired molar flows: 

^(ratio of reflux) = (dA-d* A ) 2 + (sb-s* b ) 2 

+ (b A -b* A f (2) 

where d A and d* A are the molar flows (calculated and de- 
sired) of the A (light) component in the distillation; 
Sb and s* B are the molar flows (calculated and desired) of 
the B (intermediate) component in the lateral stream; b c 



Table 2 Results obtained applying both algorithms for the pre-design of Petlyuk columns 
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Solution Description of the solution type Case of study number (see Table 1) 

type 



Table 3 Results with feasible solutions obtained with the traditional pre-design method, which uses fixed fractions 



A Feasible solution with all satisfied specifications 5, 32, 51, 70, 71, 80, 84, 89 

B Feasible solution with two output purity specifications satisfied 3, 24, 27, 33, 42, 44, 47, 50, 52, 55, 57, 59, 64, 67, 73, 78, 81, 8 

C Feasible solution with one output purity specification satisfied 4, 17, 23, 30, 38, 41, 61, 83, 90 

D Infeasible solution. There is no convergence in the pre-design 1, 2, 6, 7, 8, 9 10, 11, 12, 13, 14, 15, 16, 18, 19, 20, 21, 22, 25, ; 

Stage 29, 31 , 34, 35, 36, 37,39, 40, 43, 45, 46, 48, 49, 53, 54, 56, 58, £ 



and b* c are the molar flows (calculated and desired) of 
the C (heavy) component in the bottom stream. The cal- 
culated values refer to those obtained for the Petlyuk- 
type column and the desired values are those obtained 
in the pre-design stage. 

Once the reflux ratio is obtained, the amount of en- 
ergy used (Q) is minimized using Equation 3, defined as 
the sum of the absolute values of heat required by the 
boiler (Q r ) and the heat removed in the condenser (Q c ); 
the minimization is achieved by the variation of the 
interconnection flows (L sup , V n ') of the Petlyuk-type 
column. 

g 3 (L sa »,V M )=\Q r \ + \Q c \ = Q (3) 

The optimal design of a Petlyuk-type distillation col- 
umn is achieved after carrying out this step. 

Experimental 

The cases studied here allow us to show the benefits and 
problems of both algorithms when using them for pre- 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

F (mole frac) 

Figure 5 Composition profiles obtained for the main column 
by the traditional pre-design method. Study case 27: n-pentane/ 
n-hexane/n-heptane, Fa = 0.45, Fb = 0.10, Fc= 0.45 and q = 0.0. 



6, 87 

16, 28, 
0, 62, 

63, 65, 66, 68, 69, 72, 74, 75, 76, 77, 79, 82, 85, 88 



designing Petlyuk-type columns, and also for estimating 
general response surfaces. In all cases, it is taken for 
granted that the system is working at a constant pres- 
sure with a feed flow rate of 45.36 kgmol/h. 

In order to show the benefits of the proposed method, 
only three factors were studied. However, further studies 
should be done including more factors. The factors stud- 
ied here are: 

a) Mixture type (M). Three hydrocarbon ternary mixtures 
were studied. Two of which have been already 
mentioned in the literature regarding their separation 
in Petlyuk-type columns [3,4,17]: n-pentane/n-hexane/ 
n-heptane {Mi), n-hexane/n-heptane/n-octane (M 2 ); 

a third mixture has not been reported yet in the 
literature: n-hexane/benzene/toluene (M 3 ). 

b) Composition of the feed (Fa, Fb y Fc). Where Fa is 
the molar fraction of the most volatile component; 
Fb is the molar fraction of the intermediate volatility 
component, and Fc is the molar fraction of the smallest 
volatile component. These three quantitative factors 
represent the composition of a mixture and their sum 
must be one: Fa + Fb + Fc = 1. The composition of each 
feed was studied within a range of 0.10 to 0.80. 

c) The thermal conditions of the feed (q). This 
quantitative factor was studied within a range of 0.00 
(saturated vapour) to 1.00 (saturated liquid), that is to 
say, the mixtures were fed in a liquid-gas state. 

In the distillation columns, higher fixed costs are gen- 
erated by the column size, and higher operating costs 
are generated by thermal loads. Here we studied the ef- 
fect caused by the factors mentioned above on the fol- 
lowing two response variables: 

Table 4 Statistics of number of theoretical stages (AT) 
obtained by proposed method, grouped by mixture 



Mixture Min. Mean Max. 

M-, 40 51.77 67 

M 2 46 61.30 76 

M 3 94 130.90 159 

Total 40 81.32 159 
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Table 5 Statistics of energy consumption (Q in MW) 
obtained by proposed method, grouped by mixture 


Mixture 


Min. 


Mean 


Max. 


M, 


2.09 


7.57 


30.66 


M 2 


2.26 


9.87 


50.37 


M 3 


6.39 


56.43 


636.37 


Total 


2.09 


24.62 


636.37 



a) The total number of theoretical stages for the 
equivalent configuration (A/); this number is 
obtained as a result of an optimization carried out in 
the pre-design stage, Equation 1. The minimum 
value of this parameter was searched for; the num- 
ber of theoretical stages accounts for a large percent- 
age of the cost of equipment investment. 

b) Thermal load of the Petlyuk column (Q in MW), 
which results from solving the optimization problem 
proposed in Equation 3. A minimum value of this 
parameter was searched for, since it accounts for a 
large percentage of the equipment operating cost. 

An experimental strategy involving ninety study cases 
was established for the purpose of studying the effect that 
each of the factors mentioned may have on the response 
variables (TV and Q). The experimental plan, obtained by 
the D-optimal strategy, contains the necessary number of 
experimental points required to obtain a representative re- 
sponse surface from the combinations of the levels of all 
factors. This experimental plan ensures that we fully study 
each factor within its range, not making the mistake of 
studying only a few points that may yield good results. 
The experimental plan is presented in Table 1. An Aspen 
Plus V7.1* simulator was used in the ninety study cases for 
obtaining the necessary information in the pre-design 
stage, as well as in the optimal design stage. In all cases, 
we used the RK- Aspen property method for predicting all 
phases. The Redlich-Kwong-Aspen equation-of-state is 
the basis for the RK- ASPEN property method and can be 
consulted at [18]. 

Results and discussion 

Table 2 presents the results obtained for the number of 
theoretical stages and the thermal load when both algo- 
rithms were applied in the design of Petlyuk columns. 
Feasible results were obtained with the proposed method 
(Method 2) for all 90 study cases with a solution that met 
all the established specifications, whereas the traditional 



method (Method 1) provided feasible results only for 8 of 
the 90 study cases. 

An analysis of the failures of Method 1 is presented. 
Table 3 presents a summary of the results obtained with 
Method 1. Only 8 study cases (which represent 8.9% of 
the study cases involved) obtained feasible solutions that 
met all the established specifications of the problem 
addressed (type A solution), 5 of which are for the Mj 
mixture, 3 for the M 2 mixture, and none for the M 3 mix- 
ture. This explains, perhaps, why in previous works with 
Petlyuk columns only ideal mixtures have been used. 

Other feasible solutions were obtained for 28 study 
cases (type B and C solution), but the results were differ- 
ent from those required (type 2 solutions). For example, 
the type B solution implies that purities lower than 98% 
are obtained in a component by applying the design 
method, even when a value of 98% had been established. 
The type C solution involves a failure in the purity of 
two components. Figure 5 shows this for study case 
number 27 (type B solution), in which a purity of 55% 
was obtained for the intermediate component (n-hex- 
ane). Infeasible solutions were obtained for 54 study 
cases (60% of cases), that is to say, the FUG equations 
did not provide a solution for the conditions specified in 
the traditional method. 

In all cases where type A, B and C solutions were ob- 
tained, the number of theoretical stages in the traditional 
method was less that in the method proposed here. Al- 
though this is an advantage of the traditional method, 
the proposed method always converges to the desired 
purity specifications. Comparing these two aspects, the 
proposed method turns out to be the better one, since 
the design of the Petlyuk columns is based on the purity 
of the products. With respect to the thermal load, simi- 
lar results were obtained with both methods; in some 
cases the proposed method yielded better results than 
the traditional method and in other cases the reverse oc- 
curred. These results show that the proposed method is 
comparable to the traditional method in terms of the 
economic benefits reported in other studies, but it pro- 
vides feasible results always. Thus, the proposed method 
allows for exploring a number of conditions, which was 
not possible with the traditional method. 

Tables 4 and 5 present the descriptive statistics of the re- 
sults obtained through the proposed method (Method 2), 
regarding the number of stages and the thermal load, 
grouped by each of the mixtures. In these tables, it can be 
seen that the Mj mixture is easier to separate and that the 



Table 6 Statistics by response (number of stages and thermal load) obtained while applying an ANOVA 

Variable Transf. SSE SSR df F-value p-value R 2 ftl dJ 



N Sqrt. 3.99 335.31 23 241.10 <0.0001 0.9882 0.9841 

Q Inv. Sqrt. 0.49 1.48 23 8.68 <0.0001 0.7515 0.6649 
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M 3 mixture is the hardest to separate, since the My mix- 
ture has lower values of A? as Q and the M 3 mixture has 
much higher values. 

Since the proposed method (Method 2) was adequate in 
all studied cases, an analysis of variance (ANOVA) and a 
response surface analysis were performed using a quad- 
ratic model for the mixed factors and a linear model for 
the process factors. In this case, we used an ANOVA, 
since the response, obtained from simulations, had a cer- 
tain degree of bias, as it came from optimization results 
(Equations 1, 2 and 3) that are sensitive to initial condi- 
tions. For example, N depends on the composition of the 
flows of the interconnection streams that were obtained in 
the pre-design stage, established by optimization. This 
ANOVA is presented only to explore the response surface 
of variables N and Q. Table 6 presents a summary of the 
ANOVA performed with transformations for the number 
of stages and the thermal load. In both cases, the model is 
statistically significant (p < 0.05) with respect to the re- 
sponse and represents a good adjustment, since the coeffi- 
cient of determination (R 2 ) is in reasonable agreement 
with the adjusted coefficient of determination (R a dj)'> there- 
fore, these models can be used to navigate the design 
space. This provides us with a safe way of representing the 
response surfaces with respect to the number of stages 
and the thermal load. 



It can be observed in Figure 6 that when the light 
component (component A) is found in greater quantity, 
there is a larger number of theoretical stages and there- 
fore a higher cost in equipment. However, a lower num- 
ber of stages, and hence a lower cost in equipment, is 
achieved when the intermediate component is found in 
a larger quantity only with the Mi and M 2 mixtures 
(subfigures a-f), while the same happens with the M 3 
mixture when there is a larger quantity of the heavy 
component (subfigures g-i). With the same mixture, 
when the thermal conditions of the feed change, the 
qualitative form and the magnitude order of the re- 
sponse surface do not suffer any considerable changes. 
In this way, the cost of the Petlyuk column is little sensi- 
tive to changes in the thermal conditions of the feed. For 
the same thermal conditions of the feed, when the mix- 
ture changes the response surface undergoes consider- 
able changes, qualitative as well as quantitative. This 
implies that the capital cost of a Petlyuk column is sensi- 
tive to the type of mixture that one is trying to separate. 
Mixture M 3 is more difficult to separate (subfigures g-i). 

Figure 7 shows that a greater thermal load, and there- 
fore a higher operation cost, occurs when the heavy 
component is found in greater proportion, while a 
smaller thermal load occurs when the lighter or inter- 
mediate component is found in greater proportion. In all 




Figure 6 Response surface for the N variable as a function of the composition of the feed. Each subfigure represents a combination of 
the type of mixture with the thermal conditions of the feed. The rows represent the mixture (M h M 2 , M 3 ) and the columns represent the feed 
thermal conditions: saturated vapour (subfigures a, d and g), vapour-liquid mixture with q = 0.5 (subfigures b, e and h) and saturated liquid 
(subfigures c, f and i). 

k - 
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Figure 7 Response surface for the Q (MW) variable as a function of the feed composition. Each subfigure represents a combination of the 
type of mixture with the thermal conditions of the feed. The rows represent the mixture (M,, M 2 , Mj) and the columns represent the feed 
thermal conditions: saturated vapour (subfigures a, d and g), vapour-liquid mixture with q = 05 (subfigures b, e and h) and saturated liquid 
(subfigures c, f and i). 



mixtures, when the thermal conditions of the feed change, 
the qualitative form of the response surface does not suffer 
considerable changes. However, when the thermal condi- 
tion of the feed changes, the order of magnitude of the 
5response suffers considerable changes too, which indicates 
that the operating cost of a Petlyuk column is very sensitive 
to changes in this factor. For the same thermal conditions 
of the feed, the response surface undergoes considerable 
quantitative changes when the mixture changes. That is 
why the cost of operation of a Petlyuk column is very sensi- 
tive to the different mixtures that one is trying to separate. 

Conclusions 

The proposed strategy allow us not only to find a feasible design 
that meets output specifications but also one in which energy 
consumption is comparable to that obtained by the traditional 
pre-design model, which uses fixed fractions. We are certain 
that with this method it is possible to obtain feasible solutions 
for a wide range of conditions, allowing for more general studies 
on the design of distillation columns of the Petlyuk-type. We 
are contemplating conducting a larger study taking into account 
other process parameters (e.g. the reflux ratio). 

Nomenclature 

df Degree of freedom 



Fa Molar fraction of the intermediate volatile compo- 
nent kgmol A/kgmol 

Fb Molar fraction of the most volatile component kgmol 
B/kgmol 

Fc Molar fraction of the least volatile component kgmol 
C/kgmol 
Y-value Fisher value 
M Type of mixture 

Mi n-pentane/n-hexane/n-heptane mixture 

M 2 n-hexane/n-heptane/n-octane mixture 

M 3 n-hexane/benzene/toluene 

N Number of theoretical stages 

p-value Probability value 

q Thermal condition of the feed 

Q Thermal load MW 

Q c Heat removed in the condenser MW 

Q r heat required in the reboiler MW 

R Coefficient of determination 

.Radj Adjusted coefficient of determination 

SD Standard deviation 

SSE Sum square error 

SSR Sum square regression 
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